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ABSTRACT 


The stability of an incompressible two-fluid wheel flow to infinitesimal helical dis- 
turbances is studied. The inner fluid is heavy and has infinite electrical conductivity, 
while the outer fluid is light and is nonconducting. An axial magnetic field is externally 
imposed on both fluids. This configuration may be viewed as a Rayleigh- Taylor problem 
in the frame of the rotating fluid and is dynamically unstable. Growth rates increase 
with increasing axial wavelength and azimuthal mode numbers, but decrease with in- 
creasing axial magnetic field. By increasing the magnetic field sufficiently, the system 
can be made stable to short axial wavelength disturbances for any azimuthal mode. 



STABILITY OF TWO-FLUID WHEEL FLOWS WITH AN IMPOSED 
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SUMMARY 

The stability of an incompressible two-fluid wheel flow to infinitesimal helical dis- 
turbances is studied. The inner fluid is heavy and has infinite electrical conductivity, 
while the outer fluid is light and is nonconducting. An axial magnetic field is externally 
imposed on both fluids. This configuration may be viewed as a Rayleigh- Taylor prob- 
lem in the frame of the rotating fluid and is dynamically unstable. Growth rates in- 
crease with increasing axial wavelength and azimuthal mode numbers, but decrease 
with increasing axial magnetic field. By increasing the magnetic field sufficiently, the 
system can be made stable to short axial wavelength disturbances for any azimuthal 
mode. 


INTRODUCTION 

Vortex containment of a heavy fissioning gas surrounded by a lighter coolant gas has 
been suggested for possible application to gaseous core nuclear rockets (refs. 1 to 3). 
Another application of the vortex containment principle is found in the vortex magneto- 
hydrodynamic (MHD) generator (ref. 4). The problem investigated in this report was 
originally conceived to study the stability characteristics of the two-fluid wheel-flow re- 
actor concept of Evvard (ref. 3). In this concept, a core of heavy fissioning gas in 
solid-body rotation is surrounded by a lighter coolant gas also in solid-body rotation at 
the same angular frequency. Since the centrifugal force is directed radially outward 
from the heavy inner fluid towards the lighter outer fluid, this flow configuration is dy- 
namically unstable. When viewed from the rotating reference frame, the situation is 
analogous to the well-known Rayleigh- Taylor instability problem where a heavy fluid is 
supported against gravity by a lighter fluid. Reshotko and Monnin (ref. 5) investigated 
the nature of this two-fluid wheel-flow instability from hydrodynamic considerations only. 



They found the flow always unstable with growth rates that increase with increasing 
axial, as well as azimuthal, wave number. 

Because the gaseous core is highly ionized and therefore electrically conducting, 
Evvard suggested that an externally imposed axial magnetic field might be used to sta- 
bilize the flow. Obviously, any realistic evaluation of the stability problems involved in 
this complex flow system would have to be determined by experiment. But the proposed 
concept does suggest an idealized flow configuration whose stability can be investigated 
mathematically by a straightforward perturbation analysis. 

This idealized flow problem, treated in the present report, consists of two inviscid 
immiscible fluids of constant and uniform density separated by a cylindrical interface. 
The inner fluid is a perfect electrical conductor, while the outer fluid has the electrical 
properties of a vacuum. There is an externally imposed axial magnetic field which, at 
equilibrium, is uniform and has the same value in both fluids. This is a hydromagnetic 
flow in which the electromagnetic fields are coupled to the hydrodynamic flow fields. 

The hydromagnetic equations, with scalar pressure, were used to describe the flow. 

Although the correspondence between the real flow problem and the idealized one is 
tenuous indeed, the idealized problem possesses the advantage that it can be readily 
solved and some of the results should be of use in the design of an experiment and in the 
interpretation of results. The idealized problem is also of interest in its own right 
since it has not been previously investigated. Furthermore, in hydromagnetic fluids 
such as a liquid metal, the sharp boundary and the infinite electrical conductivity would 
correspond to a realistic physical situation. 

Wilhelm (ref. 6) has investigated the stability of an idealized two-fluid wheel-flow 
configuration where the inner fluid was a perfect conductor and the outer fluid was a 
nonconducting gas. The equilibrium magnetic field was an azimuthal field produced by 
an axial sheet current at the interface between the inner and outer fluids. The basic 
equations used in the present report are identical to those of Wilhelm. 

In the idealized model, viscous effects were omitted, but as Chandrasekhar (ref. 7) 
notes for rotating flows, viscosity severely diminishes the growth rate of short- 
wavelength disturbances. We should therefore expect the short- wavelength results of 
this report to be pessimistic in the predictions of stability. Finite resistivity effects 
have also been omitted in the analysis, but Kruskal and Schwarzschild (ref. 8) have 
found that for the kind of problem under consideration, namely large electrical conduc- 
tivity, the infinite conductivity assumption preserves the essential features of the real 
situation. 

In this report, the hydromagnetic stability problem is solved by the normal- mode 
technique. We begin by writing the continuity, momentum, Ohm’s law, and Maxwell 
equations for the inner and outer fluids. These equations are linearized by performing 
a perturbation analysis about an assumed equilibrium flow. Solutions to these equations 
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are obtained for the inner and outer fluids, and the interfacial boundary conditions are 
derived. Simultaneous solution of the volume and boundary equations yields the disper- 
sion relation. The dispersion relation is solved for complex frequencies and real wave 
numbers. The results suggest some methods that might be employed to stabilize this 
type of flow. Dimensional analysis is applied to determine the degree to which the 
idealized model is approximated by the wheel-flow reactor plasma, and by the liquid 
metals sodium, potassium, and mercury. The relevance of the idealized model to situ- 
ations where finite amplitude disturbances and finite electrical conductivity exist is dis- 
cussed in a qualitative way. 


BASIC FLOW EQUATIONS 
Inner-Fluid Equations 


The hydromagnetic equations, in MKS units, are used for the inner fluid: 


P 



+ (V* • V) V* 


= J* X B* 


VP* 


(1) 


V • V* = 0 


( 2 ) 


E* + V* X B* = 0 


(3) 


V X B* = n Q J* 
V • B* = 0 


V X E 


* 


3B* 

at 


(4) 

(5) 

( 6 ) 


The asterisks indicate that the quantities are dimensional quantities. (All symbols are 
defined in appendix A. ) In these equations viscosity, electrical resistivity, space 
charge, and displacement current have all been set equal to zero. 

The electric field E* and current density J* can be eliminated from equations (1) 
and (6) by use of equations (3) and (4) to yield (ref. 6) 
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p 


( 7 ) 


av* 

+ (v* • v)v* 
at ~ 


= — (B* • V)B* - V 




V • V* = 0 


V • B* = 0 


3B* 

V X (v X B*) = 

~ ~ at* 


( 8 ) 

(9) 

( 10 ) 


Outer-Fluid Equations 

Since the outer fluid is an electrically nonconducting fluid, the electromagnetic 
equations are uncoupled from the hydrodynamic equations. Thus, for the outer fluid, 
the hydrodynamic equations are 


VP* (11) 

( 12 ) 

The electromagnetic equations in the outer fluid are Maxwell’s equations for a vacuum, 
with the displacement current neglected. Only the magnetic field equations are re- 
quired, and these are 


av* 
) — 
at 


+ (v* • v)v* 


v . V* = 0 


V X B* = 0 

(13) 

V • B* = 0 

(14) 


BOUNDARY CONDITIONS 

Next we describe the boundary conditions of the interface between the inner and 
outer fluids. In the real physical situation, there would be a continuous transition of 
fluid properties across a thin layer separating the two fluids. As is the usual case, we 
assume the transition is discontinuous, and then allow for a sheet current and a discon- 
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tinuous jump in the field quantities. The appropriate equations may be integrated across 
the interface to obtain the following boundary conditions: 




2m 0 


= 0 


(15) 


n • (V| - Vp = 0 (16) 

n - (B* - Bj) = 0 (17) 


where ii is a unit vector normal to the interface. Equations (15) to (17) are to be evalu- 
ated at the displaced interface. A detailed derivation and discussion of the boundary 
conditions is given in appendix B. 

It is also necessary to have a relation between the velocity field and the displace- 
ment of the fluid at the interface. Since all the equations in this report are given in the 
Eulerian form rather than the Lagrangian form, we cannot, in general, express the 
various field quantities in terms of a displacement vector. However, when the displace- 
ments are infinitesimal, the Lagrangian displacement can be related to the Eulerian ve- 
locity. Following reference 9, let the surface of the interface be defined by the function 
F(r*, 0, z*, t*) = 0. Since a particle lying in the surface must move tangential to the 
surface, it follows that (ref. 9) 


°£-+ (V* ■ V)F = 0 (18) 

at* ~ 


Let the interface be described by the function 


F(r*, 9, z*,t*) = r* - [£*( 0 , z*,t*) + a] 


(19) 


where £*(0, z*, t*) is the radial displacement from the equilibrium position r* = a. 

We now obtain a relation between the interface position and the velocity by substi- 
tuting equation (19) into equation (18). This yields 



+ v; 


v 0 9g* 
r* 90 


y* 

z az* 


'At interface 


= 0 


( 20 ) 
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Equation (20) will be used to determine the perturbation velocity in terms of the radial 
displacement of the interface. 


EQUILIBRIUM FLOW 

The geometry of the two-fluid wheel flow is shown in figure 1 . Cylindrical coordi- 
nates are used, with the z- coordinate directed into the paper. The heavy inner fluid is 
given the subscript 1, and the outer fluid has the subscript 2. The interface between the 

Azimuthal 



two fluids is at radius a. The outer boundary of the outer fluid is at radius R, which is 
taken to be infinite in the following analysis. Clearly, the equilibrium flow is unaffected 
by a cylindrical boundary at a finite radius R, as long as the boundary rotates with the 
angular velocity of the wheel flow. The reason for letting R be infinite is to simplify 
the dispersion relation obtained from the perturbation analysis. Since the disturbance is 
generated at the interface and decays rapidly away from the interface, the outer boundary 
has little influence on the disturbance growth rates. Reshotko and Monnin (ref. 5) ex- 
amined the effect of the outermost boundary at r = R, and concluded that "For values of 
Pj/Pg and R of interest in wheel flow reactors, the stability characteristics are es- 
sentially those of the unbounded configuration. " 

The inner and outer fluids move in solid-body rotation at the same angular fre- 
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quency S2. Thus, the equilibrium flow velocity is azimuthal and depends only on the ra- 
dial coordinate according to the relation 


V* e = fir* (21) 

In equilibrium there are no zero-order electric currents; the externally imposed mag- 
netic field is axial and uniform throughout space. Thus, in equilibrium the flow is 
purely hydrodynamic since the magnetic terms drop out of equation (7). It was shown in 
reference 5 that V|, given by equation (21), satisfies the pertinent hydrodynamic mo- 
mentum and continuity equations (eqs. (7) and (8)). The equilibrium pressure distribu- 
tion, in the inner and outer fluids, consistent with boundary condition (18), was found to 
be 


P 0 l = p c + Pi 0 < r* < a (22a) 

* 2 


*2 2 \ 

^-\ a < r* < R (22b) 

2 / 


0,2 


= PI + 


P 1 


A 2 o2(r 

— — + p 2 n - 


FORMULATION OF STABILITY PROBLEM 
Inner-Fluid Disturbance Equations 


The equations governing the infinitesimal disturbances in the inner fluid are obtained 
from a linearization of equations (7), (8), and (10) about the assumed equilibrium flow. 
Let B* = zBq + b*, V* = §r*S2 + v*, and P = Pq + p*. Then the linearized perturbation 
equations are 

Momentum (eq. (7)): 



+ n 



- 2ftvJ = 


Bj 0bj 

P 0 P 


1 

P 




(23) 


at* 


9v* 

+ 12 h 2!2v 

de 


* 

r 


^0.^1 
P G P 9z* 





(24) 
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dvZ 3v* B* 9b* , a / B*b*' 


3t* 3 0 ii Q p 3z* p 9z 


M r 


( 25 ) 


Faraday’s law for a perfect conductor (eq. (10)): 


3b: 


9b; 


3v! 


_JL + S2 — = B* 

9t 30 0 3z* 


(26) 


3b 


e , o d *A - 


3v 


1- 06 

at* 

90 

0 9z* 

3b* 

3b* 

3v* 

~^+ n 

z 

B o ! 

at* 

90 

0 9z* 


(27) 


(28) 


Continuity (eq. (8)): 


9v* v* ! 9v* 9v* 

_j: + _i + J__? + _ z = o 


3r 


r* 30 9z 


(29) 


Equations (23) to (29) can be nondimensionalized by using the following definitions: 


V* = vv. 


ref 


t* = tn" 1 


Id 

P* = p — 

M n 


B* = BB. 


ref 


co* = cofi 


v ref = aSJ 


r* = ar 


D" 

A 2 - ref 

A a o 


^o^a v ref 


z* = az 


where subscript a = 1 refers to the inner fluid, and a = 2 refers to the outer fluid. 
The nondimens ional form of equations (23) to (29) becomes 
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(23a) 


3v r 9v r o Sb r oa 

<P+ B 0 b z, 


— + — + 2v = A^B 0 B b ) 

at 30 r 1 u az r 30 u z 


9v„ 3v „ 3b 

+ — = A ? B 0 -77 - A 1 1 to + B o b z> 


at ae 


az 


az 


fir f^r B fir 

3t 39 0 3z 


ffe + f% =B Hi 

St 39 0 3z 


3b, ab, 3v_ 
— £ + — ^ = B q — £ 
at a© u az 


3v v J 3v fl 3v 

3r r r 30 3z 


(24a) 

(25a) 

(26a) 

(27a) 

(28a) 

(29a) 


Since the coefficients of equations (23a) to (29a) are independent of t, 0, and z, the 
equations allow solutions of the form q(r)exp[i(m0 + kz + art:)] where q(r) is a complex 
disturbance amplitude, k and m are the nondimens ional axial and azimuthal wave num- 
bers, respectively, and o> is a nondimens ional complex frequency. The wave numbers 
k and m are real, and m is integral. The real part of a) is the nondimens ional ro- 
tational frequency of the disturbance, while the imaginary part o>. is the nondimen- 
sional growth rate of the disturbance. Disturbances grow, are neutral, or decay ac- 
cording to whether is negative, zero, or positive, respectively. For the assumed 
form of the disturbance, equations (23a) to (29a) become 


icv r - 2Vq = AjB Q ikb r - A* ^ + B Q 

A 2 

icVg + 2v r = A 2 B o ikb 0 im (p + B 0 b z ) 


(23b) 


(24b) 
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icv z = _A i ik P 


(25b) 


icb r = B 0 ikv r 


icb 0 = B o ikv 0 


icb z = B 0 ikv z 


dv v imv fl 

-JL+JL + - + ikv z = 0 

dr r r 


where c = u> + m, and c is the wave frequency observed in the rotating frame. It is 
subsequently referred to as the relative frequency. Substitute equations (25b), (26b), 
and (28b) into equation (23b) and equation (27b) into equation (24b), and let 


Qi = ic 1 


A 2 k 2 B§ 


Then 


Vr ' 2V 0 = 


A 1 Q 1 dp 
ic dr 


Q l v 0 + 2v r " _A 1 ““ Q l p 

JL L 


Solve equations (32) and (33) for v and v fl in terms of p to obtain 


1 / 2 dp _ im p 


6 2 ic \ Q 1 dr r 


~ 1 Q2dp + 2im Qp 
Q 2 6 2 ic V dr r 1 > 
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where 


*2 < 4 

6i=l+ — 

Q? 


(36) 


Substitute equations (25b), (34), and (35) into equation (29b) to obtain 


di + I d P_/m? + k 2 § 2\ Q 

dr 2 \r 2 7 


(37) 


The general solution of equation (37) is a linear combination of Bessel functions of 
the first and second kinds. The form used here is 


p = *^1 J m ( ik6 i r ) + J *2 H m( ik5 l r ) 


(38) 


The second term blows up at the origin, so s/ ^ is set equal to zero. The solution in 
the inner fluid is 


P 1 = J *i j m ( ik6 i r ) 


(39) 


The boxed equations will be used in the interface boundary condition equations to obtain 
the dispersion relation. We obtain by substituting p^ for p into equation (35). 


lr 


- srf A 2 

QjSjic 




2im 


J m( ik6 


ik6 1 r)| 


(40) 


The magnetic field component b^ r is obtained by substituting equation (40) into equa- 
tion (26b). 


b lr -'**1 


Q, 6,1(4 dr 


+ 2*s j m ( ik « 


,4 


(41) 


Finally, bj z is obtained by combining equations (25b), (28b), and (39). 
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( 42 ) 


Outer-Fluid Disturbance Equations 

Linearization of equations (11) to (14) about the equilibrium flow yields the disturb- 
ance equations for the outer fluid. 

Momentum (eq. (11)): 


Continuity (eq. (12)): 


Div B = 0 (eq. (13)): 


Curl B = 0 (eq. (14)): 


3v* 3v* , - * 

at* 30 0 p dr* 


(43) 


3vJ 

3t* 


■ + £2 • 


30 


+ 2S2v* = - _L He! 
r pr* 30 


(44) 


3v, 


3v* 


3 + ^ = -I3p: 

St* 30 p 3z* 


(45) 


9v* v* , dv* dv* 

— + — + ?- + — = 0 


3r 


r* 30 3z 51 


(46) 


3 b* b* , 3b* 3b* 

_£ + ^ + J__? + _z = 0 

'* r* r* 30 3z* 


3r 


(47) 


1_K 

r* 30 


3b| 

3z* 


= 0 


(48) 
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(49) 


3r* 


d K_ 

3z* 


= 0 


9b 0 + ^0 _ _1_ 1ZL = o 


3b: 


Sr 51 


30 


The nondimens ional forms of equations (43) to (50) are 


3v 3v 


at 


30 


- 2v« = 


-A 2 9p 
A 2 T" 
3r 


3v„ 3v, 


3t 


30 


+ 2v = - _£ ^ 


2 3p 

r 30 


3v„ 3V,, 0 ~ 

— £ + — 5. = -a 2 

3t 30 ^ 3z 


3v. 

3r 


j 3v fl 3v 

^ + — + i — 2 + — 5 = 0 


r 


r 30 


3z 


!5r + ?r 

3r r 


1 9b 0 9b z 

i — S + — ® = 0 


r 30 


0Z 


1 ^z 

r 30 


3b, 


3z 


= 0 


3b 3b 


3r 


3z 


= 0 


9b 0 [ b 0 


3r 


, 3b 

1 r 

r 30 


= 0 


As before, the equations allow solutions of the form q(r)exp[i(m0 + kz + cot)]. 
Equations (51a), (52a), and (53a) become, respectively, 


(50) 

(51a) 
(52a) 
(53a) 
(54a) 
(55a) 
(56a) 
(57 a) 
(58a) 
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(51b) 




Q 2 v e+ 2 v r -A|*M 


2,-i 


Q 2 v z = " A 2 ikp 


where Q 2 = ic. And equation (54a) becomes 


dv_ 

+ — + — v Q + ikv z - 0 
dr r r 


(52b) 

(53b) 


(54b) 


Solve equations (51b) and (52b) for v r and v g in terms of p to obtain 


v 0 = 


1C6 2^2 


Cdp Q 2 impN 

(dr r / 


v r = 


2_ Q dp + 2imp 

. . 2 n \ z dr r j 
1C6 2 Q 2 


(59) 


(60) 


where § 2 = * + 4 /q|. 

Substitute v 0 from equation (59), from equation (60), and v z from equation 
(53b) into equation (54b) to obtain the differential equation for the perturbation pressure 
in the outer fluid 


& + i ® - (is! + k 2 6 

dr 2 rdr \r 2 

This differential equation is of the same form as equation (37). The solutions to equa- 
tion (62) are therefore 



p = • fl/ l'W ik5 2 r ) + • j/ 2 H m( ik6 2 r ) 


(63) 
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where H. is the Hankel function of the first kind. It is usually denoted by but 

the superscript was removed to simplify notation in this report. As r — °°, the first 
term in equation (63) blows up, so is set equal to zero. The solution for the per- 
turbation pressure term in the outer fluid is therefore 


P 2 = ^2 H m( ik6 2 r) 


( 64 ) 


The perturbation velocity term v p is obtained by substituting pg for p in equa- 
tion (60), then 



(65) 


We obtain b and b from equations (55a), (56a), and (57a) which become, re- 
r z 

spectively, after assuming all first-order quantities vary as q(r)exp[i(m0 + kz + u>t)] 

db b . 

— - + — + — b„ + ikb = 0 (55b) 

dr r r U z 


im 

r 


b 


z 


ikb Q = 0 


(56b) 


^z 

dr 


ikb r = 0 


(57 b) 


Substituting b^ from equation (56b) and b r from equation (57b) into equation (55b) 
yields the differential equation for b in the outer region 



1 

r 




( 66 ) 


Again, the general solution is a linear combination of Bessel functions of the first and 
second kinds 
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(67) 


» z = *l J m< lkr > + a 2 H m< ikr ) 

As r-» the perturbation field must remain finite, so 3^ is set equal to zero. 
Then 


2z 


« 2 H m< ikr ) 


( 68 ) 


Substituting b from equation (68) into equation (57b) yields 
z 




j_ a 

ik dr 



(69) 


Application of Boundary Conditions at Interface 

The boundary conditions requiring all perturbation fields to be finite at r = 0 and 
r = °° have already been applied. The remaining conditions, applicable at the interface, 
are given by equations (15) to (17). The surface normal can be written as 


n = r + n 


The zero-order surface normal is the unit vector in the radial direction r. The 
quantity n is a first-order perturbation of the unit vector. The first-order boundary 
conditions on the normal components of velocity and magnetic field at the interface are 
obtained from equations (16) and (17) as 

r • ( v 2 - Vj) = v 2r - v lr = 0 at r = 1 (70) 

f • (b 2 - bj) = b 2r - b lr = 0 at r = 1 (71) 

Strictly speaking, these two boundary conditions should be evaluated at the displaced in- 
terface. But, to within first order, equations (16) and (17) may be evaluated at the 
equilibrium interface position, r* = a or equivalently r = 1. 

Substituting equations (35), (39), and (65) into equation (70) gives the condition of 
continuity of radial velocity at r = 1, 


16 





2 2 

|?1 J m (ik6 l) + 2imQ l J m( ik6 l>] *1 = [ Q 2 H m (ik6 2 ) + 2imQ 2 *2 (72) 

Ql 6 l ^2 6 2 

where prime denotes differentiation with respect to r. 

When equations (26b), (35), (39), and (69) are substituted in equation (71), the condi- 
tion for continuity of radial magnetic field at r = 1 becomes 

*i K J m< lk6 l> ♦ 2imQ,J m (ik6,)J = - J H^ik) (73) 

c 

The boundary condition in equation (15) states that the normal force is continuous 
across the displaced interface. Writing out equation (15) to all orders results in the fol- 
lowing equation, in physical variables: 




2M r 


2 Mr 


(74) 


where all quantities in equation (74) are evaluated at r* = a + £*. Now 


B* • B* = (B* + b* ) • (B* + b*) = B* 2 + 2Bib* + b 

— a — a \ — 0 —az) \ — 0 — az) 0 0 az 


*2 

az 


We may now write out all terms in equation (74), which yields 


Pq + P X — + P 2 ^ 2 ( a + --Ji + p*(g*) - Pj + Pl f2 

2 2 ^ 


2 (a + g*) 2 


+ Pj(5*) 


n *2 p ^ 

B o . B 0 u * 


L>o 


+ — b 2z<**> + 


b £<«*>" 


Mr 


2 Mr 


B o 2 . B o 


o J L 2 /J 'o 


+ — bj z (^*) + 




2 m. 


= 0 (75) 


The radial displacement is obtained from equation (20) where 


V* = v* 
r r 
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v£ = r*ft +V ; 


V* = v* 
z z 


Substituting these into equation (20) and treating i; * as first order yields to first order 


V 9t* r 90/ r *=a 


Assuming that goes as exp[i(m0 + o>*t* + k*z*)], we get from the last equation 


And, after rearranging, 


-io>*(; (a) + v*(a) - imO£* = 0 


S*(a) = 


v*(a) 


i(o>* + mS2) 


or 


£*(a) _ 


v» 


iS2a + m^ 


= -l 


v r (l) 


Thus, in dimensionless form the displacement is 


£(D = 


ivjl) 


(76) 


Retaining only first-order terms in equation (75) leaves 


B* B* 

p 2 Ah* + p*(|*) - Pl a 2 ^* - P *(£*) + — bj z u*) - — bf z (€*) = 0 


Mr 


Mr 


( 77 ) 


Nondimens ionalizing equation (77) gives 
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( 78 ) 


|(1) + Aj Pl (l) + A2 B() b lz (l) = ^ 1(1) + ^ a|p 2 (1) + ^ A2 B() b 2z (l) 

Pi p l p l 

o o 9 

It follows from the definition of A ^ that Aj = (p 2 /pj) Ag. K we choose V re j so that 
A^ = 1, equation (78) becomes 


m + Pi(l) + B 0 b lz d) = -* 5(1) + p 2 (l) + B 0 b 2z (l) 

P 1 


179) 


Substituting into equation (79) the quantities 5(1) from equation (76), p^(l) from 
equation (39), bj z (l) from equation (42), p 2 (l) from equation (64), and b 2z (l) from equa- 
tion (68) results in the continuity of normal force boundary condition 


J*. 


2 0 2„2 

c QjSj 


[Pl + 2imQ 1 J m (ik6 l )] + •*j J m ( lk 5i) 


- = *2 B 0 H m< ik > + •‘ / 2 H m< lk5 2> 


+ TTi [^ H n> (ik5 2> + 2imQ 2 H m< ik6 2>] < 80 > 

C Q 2 6 2 


DISPERSION RELATION 

From equation (72) for continuity of radial velocity. 





(81) 


19 



where 


* m - Q| H m< ik5 2> + 2im( W ilt6 2> 
S m = Q? Jm<« k «l) + 2imQ 1 J m (lka 1 ) 

From equation (7 3) for continuity of radial magnetic field, 

„ . s m 

2 c 2 q \ 6 \ 

Substitute equations (81) and (84) into equation (80) to obtain 


_^i_ 

^2 Q 2 6 2/ 


i-^xn 



A 


m 


2~ 2*2 
c QjSj 


k 2 B 2 Q k2B o4i H m (ik) 

W - -f V 15 ,) + ■° ' ” 

c c 2 Q2 5 fH^(ik) 


= w + 


C 2 Q 


Using the definitions of Q fl) < / m , and and rearranging terms give 




^ - 1 + k 2 B 2 
*>1 




where, after setting Aj = 1, 



( 87 ) 


D _ J m (ik5 l> + 2im 

J m ( ik6 i) Qi 


Gh m 


H’ (ik) 


H m (ik) 


q 2 sic 


60 = 1 _ — 
c 2 


( 88 ) 

(89) 

(90) 


p s H m^ ikg 2^ + 2im 

V ik6 2> Q 2 


(91) 


Qi =ic 1- 


k 2 B 2 


(92) 


5?.l 


c 2 [l-“° 


(93) 


c s u) + m 


(94) 


Prime denotes differentiation with respect to r. Thus, 




and 




r=l 
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Equation (86), with the definitions following it, is the dispersion relation. It is a com- 
plex equation since o>, and therefore c, is complex. In fact, the Bessel functions have 
complex arguments because of the 6 ff factors. The dispersion relation was solved nu- 
merically by using an iterative technique on an IBM 7094 computer. 

The dispersion relation has multiple roots of which some are nongrowing Alfven 
waves. A technique was incorporated into the program to eliminate these uninteresting 
solutions. 


APPROXIMATE SOLUTIONS TO DISPERSION RELATION 
Long Wavelength Limit 

The long wavelength solutions k — 0 to the dispersion relation have been deter- 
mined in appendix C. The results are 


Re(c) 


_ 1 _ 

9t 


or 


Re(o>) = -m + 


St 


m > 0 


and 


Im(c) - - V#m - 1 


m >0 


(95) 


(96) 


(97) 


where 


St = 


Pi 

p 2 


+ 1 


Ok.i 


p 2 


(98) 


These solutions were used as initial points in the iterative computer solutions. 
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Solutions Where Wave Frequency c = 0 


The values of k and Bq for which c goes to zero are determined in appendix D. 
Having selected a value of Bq or k, the following relation gives the value of k or Bq 
for which c = 0: 

k 2 B 2 

f = k 2 l m-l< k > K m-l« + m [Vl< k > K m< k > ‘ V«Vl»] ' m 2 l m« K m« <»»> 

l-£? 

Pi 

Equation (99) is further reduced to the limiting cases of k — 0 and k — <*> in appen- 
dix D. Equation (99) was used as a check of the computer solutions. 


DISCUSSION 

Although many stable modes of oscillation were found for the dispersion relation, 
only those solutions which indicate disturbances that grow with time are discussed. 

After the discussion of solutions to the dispersion equation, we explore some of the 
limitations of the normal mode approach to the stability problem and suggest some ways 
to interpret the results. Next, some physical implications about the magnetically sta- 
bilized plasma are considered. Potential methods are suggested to avoid the instabil- 
ities predicted by the present analysis. 

. Presented in figures 2 and 3 are curves of the imaginary and real parts, respec- 
tively, of the relative wave frequency c as a function of the axial wave number k. hi 
these plots, the azimuthal wave number is a parameter, and results are presented for 
four different values of the externally imposed magnetic field. These curves are solu- 
tions to equation (86), the dispersion relation. 

Typical fluid properties for the wheel-flow reactor (ref. 10) were used to establish 
the range of dimensionless magnetic fields found in figures 2 and 3. An angular velocity 
of 200 radians per second was selected to provide a plasma Mach number of less than 
0. 1 at the interface. A density ratio p^/pg of about 100, the largest envisioned, was 
chosen for the calculations because it is the most unstable configuration. The tempera- 
ture chosen was 40 000 K, and the pressure was of the order of 100 atmospheres. The 
reactor radius was 1 meter, and the density p^ was 10. 3 kilograms per cubic meter. 

As shown in figure 2 in all but the weakest magnetic field case, the growth rate c^ 
decreases to zero at a finite value of the axial wave number for a given value of m. 

For a given m and a given magnetic field, both the real (fig. 3) and imaginary (fig. 2) 
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Figure 3. - Disturbance relative wave frequency as function of axial wave number with azimuthal wave number as 
parameter. 


parts of c must go to zero at the same value of k. Because of the scales used in the 
figures, it appears on some plots that Cj and c rg go to zero at different values of k, 
but the computer results verify that both the real and imaginary parts go to zero at the 
same value of k. The value of k at which c goes to zero is referred to as the criti- 
cal wave number k cr ^ical' Be y° nd ^critical’ there are no solutions that correspond 
to a continuation of the curves presented in figures 2 and 3. 

In the hydrodynamic case (ref. 5), the flow was found to be always unstable. How- 
ever, figure 2 reveals that the effect of the magnetic field is to make the system stable 
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to short wavelength disturbances, A < 27r A cr itical* s ^ or ^ wavelengths are stabi- 

lized because the magnetic field lines, which are ’’frozen" into the fluid, undergo a dis- 
tortion or "stretching" when the interface is perturbed. It requires energy to stretch 
the field lines. Thus, the energy that would normally be available to drive the instabil- 
ity is absorbed in the stretching process. In the hydrodynamic Rayleigh- Taylor prob- 
lem, surface tension provides the "stretching” mechanism that stabilizes the short- 
wavelength disturbances. 

Figure 2 also reveals that for a given value of k, the growth rate increases as the 
azimuthal wave number increases. This suggests that smaller chunks of the heavy fluid 
have the most likelihood of tearing off at the interface, hi figure 2, the m = 1 curves 
have a dip in them near k = 0. Although it is not obvious from the figure, the growth 
rate is finite at the minimum of the dip. This same peculiar dip occurred in refer- 
ence 5, and we have no good physical explanation for its occurrence. 

We note from figure 2 that the k = 0 disturbances (infinitely long wavelength) are 
unstable for all azimuthal wave numbers except m = 0. (When m = 0 and k = 0, there 
is no disturbance on the flow. ) These unstable modes are the pure flute modes. They 
do not distort the magnetic field and, hence, are not stabilized by it. When k = 0, the 
dispersion relation reduces to the simple form given in equations (95) to (98). From 
equation (95) and from figure 3, note that c re reduces to 1 /St , when k = 0. Thus, to 
an observer moving with the rotating frame, all flute modes oscillate with the same 
frequency. 

In figure 4, the value of magnetic field for which c goes to zero is presented as a 
function of k cr y_ ical , with m as a parameter. These are plots of equation (99) in 
which the dispersion equation was solved subject to the requirement c = 0. These 
curves have a simple interpretation: for those combinations of magnetic field and axial 
wave number that lie above the curve, the plasma is stable, while it is unstable if the 
combination falls below the curve. Information found from figures 2 to 4 and equa- 
tion (C9) indicates that no matter what value of magnetic field is chosen, a sufficiently 
high value of m will make the system unstable to long- wavelength disturbances. But 
when the value of m becomes so large that the azimuthal wavelength is comparable 
with the boundary layer thickness, the analysis breaks down. The analysis implies that 
even the idealized concept of a two-fluid wheel-flow device is impractical unless some 
way is found to place restrictions on the axial wavelength and azimuthal wave numbers. 

The dispersion relation was obtained from the linearized perturbation analysis about 
an assumed equilibrium flow. This suggests several interesting limitations to the anal- 
ysis. We shall address ourselves to three of these problem areas. 

First, it was assumed in the theoretical analysis that a particular equilibrium flow 
could be established. Experimentally this may not be possible. However, the stability 
analysis gives us some insight about the assumed equilibrium flow. If in establishing 
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Figure 4. - Neutral stability curves for various azimuthal wave numbers. (Flow is stable above curve and unstable 
below. ) 


the equilibrium state the flow must pass through some of the perturbations that are pre- 
dicted to be unstable, we can be certain of great difficulties, if not the impossibility, in 
reaching the equilibrium situation. However, the analysis suggests that if the disturb- 
ances could be restricted to short axial wavelengths, the equilibrium flow might be es- 
tablished. This might be accomplished by keeping the axial length of the machine short. 

The second limitation is concerned with the assumption of small amplitude disturb- 
ances and the attendant linearized analysis. Unstable solutions of the linearized prob- 
lem represent only the first stage of growth, or of departure, of the system from the 
equilibrium state. As the disturbance amplitude becomes large, the flow should be de- 
scribed by nonlinear differential equations. It is possible for nonlinear effects to limit 
the disturbance amplitudes and produce nongrowing, finite amplitude disturbances. We 
do not expect this in the present case. The state of minimum potential energy for the 
wheel-flow problem corresponds to the heavier fluid being at the outside with the lighter 
fluid at the center. The nonlinear effects may control the growth rate of the instabili- 
ties, but should not terminate the growth. Even in the case of the short- wavelength 
disturbances, the predicted stability applies only to small amplitudes. For large am- 
plitude disturbances at these wavelengths, the nonlinear effects could produce a positive 
growth rate. 
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The linearized perturbation analysis is useful in two respects. First, of course, it 
warns of regions where the system is unstable even to the infinitesimal perturbations. 
Second, it suggests which system parameters might be manipulated to maintain or im- 
prove the stability in regions where the infinitesimal perturbations do not grow. The 
linear analysis shows that improved stabilization is obtained by increasing the magnetic 
field and restricting the disturbances to short wavelengths. Stabilization against large 
amplitude disturbances might be expected to improve with use of even larger magnetic 
fields and shorter wavelengths. The utility of this interpretation was demonstrated in 
reference 11, where the stability of a liquid- mercury - air interface was experimentally 
investigated. In these tests, mercury in glass tubes was subjected to vibrations of var- 
ious frequencies and amplitudes, and the limits of interfacial stability were noted. As 
the amplitude and frequency of vibration were increased, a smaller tube diameter was 
needed to stabilize the interface. This was consistent with the Rayleigh- Taylor theory, 
with surface tension effects included, which predicts that the interface is stable against 
short-wavelength disturbances. 

Finally, we want to know how the results of this report would be altered by the in- 
clusion of finite electrical conductivity. Kruskal and Schwarzschild (ref. 8) obtained an 
approximate solution to a simple hydromagnetic problem with finite conductivity and 
small perturbations. They compared the results with the infinite conductivity case. 

They concluded that the essential features of the large but finite conductivity problem 
were fairly represented if one assumed infinite conductivity, while allowing for electric 
sheet currents and electric sheet charges at the plasma surface. 

Intuitively, it seems that the result of finite, but large, conductivity would be some- 
what larger growth rates and less stability at the shorter wavelengths than predicted by 
the infinite conductivity case. The effect of finite conductivity is to allow the fluid to 
slip through the magnetic field, without a proportional distortion of the field. The 
energy stored in the distorted field is less than it would be for infinite conductivity. 

Finite conductivity is accompanied by a dissipative effect because a fraction of the field 
energy created by the magnetic field distortions is converted into joule heat. By con- 
trast, the energy- storage capacity of surface tension is independent of the growth rate 
because it is not accompanied by any dissipation of energy. 

The importance of the finite conductivity can be deduced from dimensional analysis. 
When the fluid has finite conductivity, equation (10) is modified by the inclusion of 
another term, becoming 


1 9 0D 

B* + V X (V* x B*) = 

" at* 


When this equation is nondimensionalized, we get 


(100) 
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1 


( 101 ) 


9 0B 

V^B + V X (V X B) = — 

fia 2 a/n o 9t 

o 

This equation can be reduced to equation (10) if S2a on » 1. The dimensionless 

2 ^ 
quantity S2a oi± Q is the magnetic Reynolds number. When the magnetic Reynolds num- 
ber is large, the field restrains the fluid as it does in the infinite conductivity case, 
where fluid stays ’’frozen” to the field lines. The criterion of large magnetic Reynolds 
number can be rewritten as a requirement on the electrical conductivity, namely 


fia 2 ^ 

-7 

Choosing O = 200 radians per second, a = 1 meter, and ju Q = 47rXlO henry per meter, 
we find that o » 4000 mho per meter is the requirement for the infinite conductivity 
assumption to apply. 

For the proposed wheel-flow reactor, the plasma has a temperature of about 
40 000 K and a pressure of a few hundred atmospheres. To estimate the electrical con- 
ductivity, we assume that the plasma is fully ionized and use the formula given by 
Spitzer (ref. 12, eq. (5-37)) for electrical conductivity. In deriving the expression for 
conductivity, Spitzer neglected terms of the order 1 /In A, where In A is the Coulomb 
logarithm. In the proposed wheel-flow reactor, In A is of the order of 2. 5 to 5. Thus, 
the conductivity can be estimated to an accuracy of 20 to 40 percent at best. The pres- 
ence of multiply ionized ions would increase the calculated conductivity by slightly less 
than a factor of 2, while partial ionization would decrease the conductivity due to 
electron-neutral collisions and recombination. But the calculated conductivity should 
be valid to within an order of magnitude, which is sufficient for our purposes. 

Based on Spitzer's formula, the conductivity is between 2500 and 5000 mho per 
meter. Hence, finite conductivity effects will play an important role in the proposed 
wheel-flow reactor. The flow could, therefore, be substantially more unstable than 
predicted by the infinite conductivity analysis. 

In contrast to the wheel-flow reactor plasma, liquid metals have much higher con- 
ductivities. Typical conductivities for the liquid metals sodium, potassium, and mer- 

7 fi fi 

cury are respectively 10, 7x10 , and 10 mho per meter. These conductivities are 
well above 4000 mho per meter, so the infinite conductivity approximation should apply 
to two-fluid wheel-flow configurations, where the inner fluid is one of these liquid 
metals. 

A second criterion for the infinite conductivity analysis to apply involves the char- 
acteristic time in which magnetic field perturbations are dissipated by joule heating. 
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When the conductivity is so low that the diffusion term in equation (100) dominates over 
the flow term, we obtain the diffusion equation 


ctm 0 « 

By dimensional arguments, the characteristic time for magnetic field inhomogeneities 
to decay or smooth out is shown to be 


* 2 
r = a a m 0 

For disturbances which occur in times short compared to r*, the fluid responds as 
though it had infinite conductivity. The growth times in this report, and in reference 5, 
are proportional to 1/c^l. For the infinite conductivity results of this report to apply, 
the growth times should be short compared to t*. Thus, the second criterion becomes 

2 

r*flc i = Oa a/i Q c. » 1 


While the first criterion required only that the magnetic Reynolds number be large, 
the second criterion requires that the product of the magnetic Reynolds number and the 
growth rate must be large. From figure 2, the growth rates Cj are of the order of 3 
or less for m < 5. In applying the second criterion, the growth rate Cj shall be given 
the value 3.0. 

For the plasma and liquid- metal systems of interest, the following table gives 
typical values for magnetic Reynolds number, magnetic decay time, and the ratio of 
magnetic decay time to growth time. 



Magnetic 

Magnetic 

Ratio of magnetic 


Reynolds 

decay time, 

decay time to 


number, 

T*. 

growth time, 


2 

fia a/i Q 

sec 


Plasma 

0.6 

3xl0~ 3 

1. 8 

Sodium 

2520 

12. 6 

7560 

Potassium 

1680 

8.4 

5040 

Mercury 

260 

1.3 

780 
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For the wheel-flow reactor, both the magnetic Reynolds number and t*£2Cj are of 
order 1, implying that the magnetic field will have very little influence and that the sta- 
bility of such a system would be better represented by the hydrodynamic analysis of ref- 
erence 5. 

For the liquid metals, the opposite conclusion applies. The magnetic effects 
clearly dominate. Magnetic stabilization of a free surface may have application to those 
space power conversion systems that employ a liquid metal as the working fluid in a 
Rankine power cycle. A stable liquid- vapor interface must be maintained in the con- 
denser to prevent vapor-locking of the condensate pump. Although there is a zero- 
gravity environment, there are local mechanical disturbances that tend to destabilize 
the interface. The results of this report, and those found in reference 13, indicate that 
strong magnetic fields and short wavelengths will provide a stabilizing influence. The 
disturbances can be restricted to short wavelengths by making the boundaries small. If 
the confining walls are made of metal, induced voltages in the fluid will be short- 
circuited. And like a short-circuited homopolar generator, the movement of the fluid 
would be severely restricted. 


SUMMARY OF RESULTS 

The stability of an incompressible two-fluid wheel flow to infinitesimal helical dis- 
turbances has been considered. The inner fluid is heavy and has infinite electrical con- 
ductivity, while the outer fluid is light and is nonconducting. Both fluids are inviscid. 

A constant and uniform axial magnetic field is impressed on both fluids. Important re- 
sults are: 

1 . The growth rates of both propagating and nonpropagating waves are diminished 
by increasing axial magnetic field. 

2. For specific values of axial and azimuthal wave number, a value of magnetic 
field can be found which reduces the growth rates to zero. 

3. Growth rates at a specific axial wave number increase for waves with increasing 
azimuthal wave number. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 23, 1969, 

129-02-08-02-22. 



APPENDIX A 


SYMBOLS 


A 

area 

L 

length 

s4 

constant in solution of differ- 

m 

azimuthal wave number 


ential equation 

n 

vector normal to interface 

a 

interface radius 


surface 

B 

magnetic field 

n 

first-order perturbation of unit 

3 

constant in solution of differ- 
ential equation 

P 

vector normal to surface 
pressure 

b 

disturbance magnetic field 

P 

disturbance pressure 

c 

wave frequency observed in 

Q 

defined by eq. (89) or (92) 

c i 

rotating frame, u> + m 
growth rate 

q(r) 

depending on subscript 
complex disturbance amplitude 

D 

defined by eq. (87) 

R 

outer radius 

E 

electric field 

Re 

real part of 

F 

defined by eq. (91) 

a 

defined by eq. (98) 

G 

defined by eq. (88) 

s 

ratio of Alfven speed to refer- 

H m (x) 

Hankel function of first kind, 


ence speed 


argument x 

s 

surface variable 

X ' 
m 

defined by eq. (82) 

t 

time 

LJx) 
m v ' 

modified Bessel function of first 

V 

velocity 

kind, argument x 

V 

disturbance velocity 

Im 

imaginary part of 

W{a, b} 

Wronskian of a and b 

J 

current density 

r, 0,z 

cylindrical coordinates 

J„(x) 

irr 

Bessel function of first kind, 

x,y, z 

Cartesian coordinates 

argument x 

6 

defined by eq. (90) or (93) de- 

m 

defined by eq. (83) 


pending on subscript 

K m< x ) 

modified Bessel function of sec- 

% 

radial displacement of inter- 

k 

ond kind, argument x 
axial wave number 


face from equilibrium 
permeability of free space 
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p 

density 

ref 

reference quantity 

O' 

electrical conductivity 

a 

dummy script 

ft 

angular wheel-flow velocity 

0 

equilibrium quantity 

Cl) 

complex angular disturbance 

1 

inner fluid 


velocity 

2 

outer fluid 

Subscripts: 


vector 

c 

equilibrium quantity at r* = 0 

Superscripts: 

critical 

value of k for which c goes 
to zero 

t 

prime denoting differentiation with 
respect to r 

i 

imaginary 

* 

dimensional quantity 

re 

real 


unit vector 
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APPENDIX B 


DERIVATION OF BOUNDARY CONDITIONS 

The boundary condition equations (eqs. (15) to (17)) are found by integrating appro- 
priate differential equations across the interface. In the real physical situation, there 
is a rapid but continuous transition of fluid properties across a thin layer separating the 
two fluids. Therefore, the differential equations that we use must be general enough to 
describe the flow on either side of, and through, the transition region. Therefore, we 
assume that both fluids satisfy the following conditions: the conductivity is finite; vis- 
cosity, space charge, and displacement current can be neglected; and the magnetic 
permeability is that of free space. Any idealized fluid properties that are different on 
opposite sides of the interface, such as infinite conductivity on one side and zero con- 
ductivity on the other, must be introduced into the boundary conditions after the limits 
are taken. 

Following the approach used in reference 14, we form a volume element, fixed in 
space, which instantaneously includes the interface (fig. 5). The volume will be a right 
cylinder such that at time t the lower face of the cylinder coincides with the interface. 
As time increases to t + 5t, the interface moves up to the position shown in figure 5. 
The length of the cylinder A is considered to be very small. Thus, in the limit, the 
side walls make no contribution to those surface integrals with finite integrands. The 
surfaces at the top and bottom of the cylinder are small enough in area A to allow a 



Figure 5. - Volume element intersecting the inter- 
face between regions 1 and 2. 
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H 


surface integral to be approximated as the product of the integrand and the surface area. 

The use of a fixed volume in space precludes the difficulties associated with moving 
reference frames. When integrating across the interface, while in a frame moving with 
the interface, the field quantities in the differential equations would have to be trans- 
formed into those fields that the moving observer would experience. The resulting 
boundary conditions for the moving frame then have to be rewritten in terms of the 
fields observed in the laboratory frame. 


Derivation of Equation (16) 

Equation (16) is obtained by integrating the conservation- of- mass equation over the 
volume element of figure 5. The integral form of the conservation of mass over a fixed 
volume is 


_9 
af'v 


f pdv + f V • (pV)dv = 0 


(Bl) 


Using Gauss' theorem to transform the second integral of equation (Bl) into a surface 
integral, we get 


_a 

at^v 


/ p dv + f pV • n ds = 0 


(B2) 


The first integral in equation (B2) is approximated as 


d_ 

at 


•/ pdv = «{/ (p W dT -/ (p) t dv } 


— f p dv = — i (p,\ V • n6tA + (p 0 \ (A - V • n6t)A - (p 9 ) AA l 

3t“V 6t[ V 1/ t+6t~ V ZJ t+5t ~ V J 


Using a Taylor expansion, we get 


(B3) 


(B4) 


H** = H + 


'V 
^ at > 


6t + . . . 


(B5) 


35 


(B6) 





Substitute for (p-^j and from equations (B5) and (B6), respectively, into 

t+ 5t t - ! - 5t 

equation (B4) to obtain 



Now passing to the limit we obtain from equation (B8) 


lim — f p dv = (p t - p 9 ) V • nA (B9) 

8t, A— 0 St “v 1 


The second integral in equation (B2) can be approximated as 


J' pV • n ds = p^Vj • (-n)A + P 2 X.2 ’ nA + ^(A) (BIO) 


where <?(A) is the contribution from the cylindrical surface. This term will vanish in 
the limit A — 0 because the integrand is finite over the cylindrical surface. Thus, the 
limit of equation (BIO) becomes 

lim f pV • n ds = (p 9 V 9 • n - p 1 V 1 • n)A (Bll) 

A— 0 “s _ ^ * L 

Substituting equations (B9) and (Bll) into equation (B2) gives 

(Pl - P 2 ) V • n + p 2 V 2 • n - p t Vj • n = 0 
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And on rearrangement, we have 


Pl(V • n - • n) - p 2 (V • n - Vg • n) = 0 


(B12) 


Equation (B12) is satisfied if 


V • n = Vj • n = V 2 • n 


(16) 


Equation (16) is the kinematic condition that the fluid velocity component normal to 
the boundary must equal the velocity of the boundary normal to itself. This applies to 
our situation, where the two fluids are immiscible. In this case no mass flows across 
the surface, and the fluids on opposite sides of the surface remain in contact with the 
surface. 

When equation (B12) is satisfied but equation (16) is not, we have the boundary con- 
dition for a shock or detonation front where fluid can cross the interface. 


Derivation of Equation (17) 

Equation (17) is obtained by integrating the V • B = 0 equation (eqs. (9) and (14)) 
over the volume element of figure 5. 


J* V • B dv = 0 


(B13) 


Using Gauss' theorem to convert to a surface integral, equation (B13) becomes 


f B • n ds = 0 


(B14) 


The surface integral can be approximated by 

f B • n ds = B . • (-ii)A + B« • iiA + ^(A) = 0 
L 

And passing to the limit gives 


(B15) 
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Derivation of Equation (15) 


Equation (7) is applicable to both sides of the interface, subject to the conditions 
noted at the beginning of this appendix. Hence, we integrate equation (7) over the volume 
element of figure 5. 



+ (V • V) V 


dv = 


— /* (B • V )Y dv " 
M 0 V 



(B17) 


By making use of the continuity equation and the identities given below, we can replace 
the volume integral on the left side of equation (B17) by one that is more convenient to 
evaluate. Accordingly, we write the differential form of the continuity equation as 


ij2 = -v . (pv) 
at 


We also have the identity 


P 


av 

at 


at ~ at 


(B18) 


(B19) 


If we substitute dp/dt from equation (B18) into equation (B19), we get 

av a(pv) 

p — = — + V • V(pV) (B20) 

at at 

Using equation (B20), the integrand of the volume integral on the left side of equa- 
tion (B17) becomes 



I 


9V a 

p — + p(V . V)V = — (pV) + p(V . V)V + V V . (pv) (B21) 

at ~ - at “ - -- - 


or 


av a 

p + p(V • V)V = — (pV) + V • (pW) 

at ~ ~ at ~ 


(B22) 


Substituting equation (B22) into equation (B17) gives 


/ — (pV)dv + / V • (pVV)dv = -±- f (B • V)B dv - / V(I 

at v M 0 y v J \ 


P + — dv (B23) 
2M r 


By application of the various forms of Gauss’ theorem, the last three integrals in equa- 
tion (B23) can be converted into surface integrals. Thus, 


f v • (pVV)dv = f pVV • n 


ds 


(B24) 


f (B . V)B dv = f BB • n 
•At 


ds 


i 


2 


/( 


V|P + ^— ]dv= / (p + J^-lhds 


(B25) 


(B26) 


'K 


Substituting equations (B24) to (B26) into equation (B23) gives 


_a_ 

at 


f (pV)dv + f pVV • n ds = — f BB • n ds - / (l 

v ^ J \ 


P + — In ds (B27) 


2 K 


Each integral in equation (B27) shall be separately evaluated for the volume element of 
figure 5. 

The first integral in equation (B27) is evaluated in the same manner as for the first 
integral in equation (B2). Thus, we can write the result by replacing p in equation (B9) 


39 



by pV, and we get 


lim A f pV dv = ( Pl V* - p 9 V 9 )V • nA 
6t, A— 0 at J y - 1-1 - 


The second integral in equation (B27) is approximated as 


f pVV • n ds = P 2 V 1 V 1 • (-n)A + pgV^ • nA + /7(A) 


Passing to the limit gives 


lim y pVV • n ds = (p 2 y_ 2 y_ 2 • n - pj Vj Vj • n)A 


We may use equation (16) to reduce equation (B30) to 


lim y pVV • n ds = (p 2 V 2 - pj Vj)V • iiA 


The third integral in equation (B27) is approximated as 


y BB • n ds = B^Bj • (-n)A + BgBg • nA + ^(A) 


Passing to the limit gives 


lim f BB • n ds = (B 9 B 9 • £ - B 1 B 1 • £)A 
A~0*4 ~ l ~ l 


By using equation (17), we can simplify equation (B33) to 


lim / BB • n ds = (B 9 - B,)B • nA 
A-O y s 1 


where 


B • n = Bj • n = B 2 • n 


(B28) 


(B29) 


(B30) 


(B31) 


(B32) 


(B33) 


(B34) 


(B35) 
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Finally, the fourth Integral in equation (B27) is approximated as 



Passing to the limit gives 



(B36) 


(B37) 


Now substitute equations (B28), (B31), (B34), and (B37) into equation (B27) to obtain 


(B38) 

The two terms on the left side of equation (B38) cancel each other, while the two terms 
on the right side are orthogonal to each other. The orthogonality condition follows from 
the fact that the vector Bg - B^ lies in the surface of the interface. This is easily seen 
by noting that the magnetic field vectors Bj and Bg may be resolved into components 
normal and tangential to the interface, B n and B^ respectively. Thus, 



(Pi Xi - p 2 V 2 ) V • n + (p 2 V 2 - pjVj) V • 5 = (B 2 - B^B • £ - n 



—2 


*1 


*n2 


5*2 


- B 


hi 


- 5*1 


(B39) 


But from equation (17), the equation of continuity of the normal component of magnetic 
field at the interface, we have that 


5 n2 


= B 


hi 


(B40) 


Thus, 


-2 " -1 “ ^t2 " ^tl 


(B41) 


Therefore, each of the two terms on the right side of equation (B38) must independently 
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be zero. This results in the following two boundary conditions: 



( 15 ) 


(B 2 - BjKB • n) = 0 


(B42) 


Equation (B42) is not used in the text because it does not yield any zero-order or first- 
order information for the problem treated in this report. However, equation (B42) can 
be interpreted in general. It is interesting to note that if B 2 - B^ ^ 0 (i. e. , when there 
are surface currents, which in turn means that one of the fluids is a perfect conductor) 
equation (B42) requires that 


B'n = Bj*n = B 2 'h = 0 (B43) 

We infer from equation (B43) that the field on both sides of the interface is wholly 
tangential to the interface. This follows physically from the infinite conductivity as- 
sumption which causes the field lines and the fluid to be "frozen” together. We know 
from equation (18) that a particle, once on the surface, always remains on the surface. 
Since the magnetic field lines and particles are always "frozen" together, we should 
expect that the magnetic field lines, once on the surface, will remain always on the sur- 
face. 

E Bq2 ^ 5oi’ equation (B42) yields first-order information about the boundary con- 
ditions. To demonstrate this we write out the zero- and first-order terms in equa- 
tion (B42). 


(B 2 - Bj)(B • n) = 0 (B44) 

(log + b 2 - Bqj - ^[(B^h- b ft ) •(£ + *)]=<> (B45) 

where a = 1,2 

(Boa + !>2 ' ®01 - * +!o a -!t + iV f + V2>* 0 I 846 ) 
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H 


The zero-order terms of this equation are 

(§02 ' 2 ol>®Oa • ?) = 0 


(B47) 


The first-order terms are 

(§02 - ®01><5o<* • « + b a • f) + Q b 2 - bjXB^ • f) = 0 (B48) 

If ®02 ^ loi* Moq, * r = 0 from equation (B47), and equation (B48) reduces to 

<5o2-5oi>®o a -2 + V f > = 0 < B49) 

and since B^ 4 Bqj, equation (B49) requires that 

B n * n + b • r = 0 a = 1, 2 (B50) 

— 0 a — —a ’ 

Equation (B50) is the first-order boundary condition that results when 
the problem considered in this report, Bgg = Bqj, and thus equation (B50) does not 
apply. Furthermore, in this report *0a't = ° ; hence, equation (B48) yields no first- 
order information. 


Other Boundary Conditions 

In reference 6, Wilhelm presents a boundary condition (eq. 1.7 of ref. 6) that was 
derived by integrating equation (10) of this report across the interface. Now equa- 
tion (10) was obtained by substituting equation (3), which holds only for infinite conduc- 
tivity, into equation (6). Since equation (10) does not hold on both sides of the interface, 
it cannot be used to obtain boundary conditions. Fortunately, equation 4. 6 of reference 6 
can be obtained from equation (B38), so Wilhelm’s results are valid. 

hi summary, the following table shows the differential equations and the corre- 
sponding boundary conditions obtained by integrating the equations across the interface. 
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APPENDIX C 


LIMITING SOLUTION OF DISPERSION RELATION FOR 
SMALL AXIAL WAVE NUMBER 


The dispersion relation equation (42) is to be solved in the limit of small axial wave 
number k. The Bessel functions become 


fc J -n ( H r= , 

v ik 6 i> 

H m (lk6 2 r ) 

fc H " (ikr) L 


= m + . . . 


= -m + . . . 


H m (ik) 


= -m + . . . 


(Cl) 


(C2) 


(C3) 


For k — 0 


*2 *2 _ « 4 

5 1 - 6 2 - 1 " 

c 4 


Qi - Q 2 = ic 


(C4) 

(C5) 


The dispersion relation can be rewritten 


(c 2 - 4)(<*c 2 - 2c + m) = 0 


(C6) 
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where 


St = 


fi + i 

9 _2 

* -i 
p 2 


The solutions to the dispersion relation are 

c = ±2 


(C7) 


(C8) 


c=lfl ± i\&TT) (C9) 

SI ' ' 

2 2 

For the solutions given in equation (C8), = 6 2 = 0. These solutions are rejected 

because the condition that the real part of k6 2 is greater than zero does not hold. 

The solutions in equation (C9) are both valid solutions. The solution with the nega- 
tive imaginary root grows while the solution with the positive imaginary root decays. 
Although both are possible, the growth solution is of more interest here since it persists. 
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APPENDIX D 


DETERMINATION OF THE POINT WHERE WAVE FREQUENCY c = 0 


Let c = 0. The dispersion relation (eq. (86)) becomes 


0 = — - 1 + k 2 B? 

Pi 


J m< ik > H m( ik > 


H m( ik ) 


L J m< ik > 


(Dl) 


Differentiating the Bessel functions and rewriting them as modified Bessel functions 
gives 


1 - — = k 2 B 2 

P 1 


I (k) 

m v ' 


+ — - 


K m< k > 


L kI m _i( k ) - ml m (k) mK m (k) - kK m _ l( k) 


(D2) 


1 - — = k 2 B 2 . 

Pi 


k W k ) K m-l( k ) + kI m-l( k ) K m< k > 


l k2l m-l( k ) K m-l( k ) + 

(D3) 


From the Wronskian of the modified Bessel functions 

W [ K m-l«> Vi®] “ Vl< k > K m< k > - W k > K m-l< k > = £ <«> 

K 

Thus, the numerator of the factor containing Bessel functions equals 1. Rearranging 
terms gives 

k 2 B 2 

+ mk [Wl< k > K m< k > ‘ U k ) K m-l«l ‘ ” 2 W k > K m< k > 

1 - 

Pi 


(D5) 

Solving equation (D5) for B z yields the value of k at which c = 0. For several limit- 
ing cases, the equation can be further simplified. For all m as k - °°, equation (D5) 
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I III III ■ 


I I ■! II INI 


reduces to 


B 


2 

0 



Pi 


2k 


(D6) 


For m = 0, k — 0, equation (D5) becomes 



P 1 


For m = 1, k — 0, equation (D5) becomes 


R 2 

_P2 2k 2 2 

P 1 


For m > 1, k — 0, equation (D5) becomes 

R 2 

**0 1 | m 

1 P2 4 ( m " 1 ) 2k 2 

Pi 


(D7) 


(D8) 


(D9) 
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